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Abstract. Several solar analogs have been identified in the library of high resolution stellar spectra taken with 
the echelle spectrograph ELODIE. A purely differential method has been used, based on the \ 2 comparison of 
a large number of G dwarf spectra to 8 spectra of the Sun, taken on the Moon and Ceres. HD 146233 keeps 
its status of closest ever solar twin (Porto de Mello & da Silva. IT997II . Some other spectroscopic analogs have 
never been studied before, while the two planet-host stars HD095128 and HD186427 are also part of the selection. 
The fundamental parameters found in the literature for these stars show a surprising dispersion, partly due to 
the uncertainties which affect them. We discuss the advantages and drawbacks of photometric and spectroscopic 
methods to search for solar analogs and conclude that they have to be used jointly to find real solar twins. 
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1. Introduction 

The Sun is the best-known star : its fundamental parame- 
ters (radius, mass, age, luminosity, effective temperature, 
chemical composition) are known with a good accuracy, 
as well as its internal structure, activity, velocity field and 
magnetic field. Consequently the Sun is used as the funda- 
mental standard in many astronomical calibrations. One 
of the motivations to identify stars that replicate the solar 
astrophysical properties is the necessity to have other ref- 
erence stars, observable during the night under the same 
conditions as any other target. The pioneers of the subject 
(Hardorp 1978 , Cayrel de Strobel et al 193J were also in- 
volved in resolving the problem of the photometric indexes 
of the Sun, inherent to the impossibility to observe it as a 
point-like source. In the last decade the motivation of find- 
ing such stars has been increased by an exciting challenge 
: the search for planetary systems that could harbour life. 
Solar analogs are straightforward targets for this hunt. 

The first searches of solar analogs were performed by 
photometric and spectrophotometric techniques. Hardorp 
l|1978f) compared UV spectral energy distributions of 
nearly 80 G dwarfs to that of the Sun and found 4 stars 
that had a UV spectrum indistinguishable from solar : 
HD028099 (Hy VB 64), HD044594, HD186427 (16 Cyg B), 
HD191854. Neckel ljTO5rJ|l established a list of bright stars 
with UBV-colours close to those of the Sun and confirmed 
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the photometric resemblance of Hy VB 64 and 16 Cyg B 
to the Sun. With the advance of techniques in high res- 
olution spectroscopy and solid state detectors, and with 
the progress in modelling stellar atmospheres, measure- 
ments of (T c ff, logg, [Fe/H]) became of higher precision 
allowing the search for solar analogs by comparing their 
atmospheric parameters to those of the Sun. G. Cayrel de 
Strobel made a huge contribution to the subject with the 
detailed analysis of many candidates (Cayrel de Strobel 
et ai rtWn Cayrel de Strobel & Bentolila 1989, Friel et 
al 1993) and a review of the status of the art (Cayrel de 
Strobel 1996|. She also introduced the concepts of solar 
twin, solar analog and solar-like star. Porto de Mello & da 
Silva l(TM7|> presented the star HD 146233 (18 Sco) with 
physical properties extremely close to those of the Sun, as 
the " closest ever solar twin" . 

A workshop on Solar Analogs was held in 1997 at the 
Lowell Observatory to provide a solid basis to the hunt 
of solar analogs. After many discussions on the perfor- 
mances of different methods, a list of the best candidates 
was established, in which 4 stars received the agreement 
of almost all participants : HD217014 (51 Peg), HD146233 
(18 Sco), HD186408 (16 Cyg A), HD186427 (16 Cyg B). 

In this paper, we take advantage of a large and homo- 
geneous dataset of high resolution echelle spectra which 
are compared directly to solar spectra, independently of 
any model or photometric measurements. The eye is re- 
placed by a more reliable criterion, approximatively the 
reduced \' \ computed over ~ 32000 resolution elements. 
This purely differential method allowed us to identify sev- 
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eral stars whose optical spectrum looks closely like the 
Sun's, the best one being HD146233. 

We describe in Sect. our observational material and 
differential method, and we give the list of our Top Ten 
solar analogs. We have searched the literature for their 
colour indexes and atmospheric parameters and calculated 
absolute magnitudes from Hipparcos parallaxes. We dis- 
cuss the uncertainties which affect these data and compare 
them to that of the Sun (Sect. |3J). In Sect. 0] we exam- 
ine qualitatively their Li content and give information on 
their activity and age. In Sect.[S]we discuss several stars, 
having similar colours and absolute magnitude or similar 
atmospheric parameters as the Sun but slightly different 
spectra. 

2. ELODIE spectra and the TGMET code 

All the spectra used in this paper were extracted from the 
library of stellar spectra collected with the cchcllc spectro- 
graph ELODIE at the Observatoire de Haute-Provence by 
Soubiran et al ijlMSjl and Prugniel & Soubiran (j!Ml"j) . The 
performances of the instrument mounted on the 193 cm 
telescope, are described in Baranne et al l|1996). ELODIE 
is a very stable instrument, built to monitor radial ve- 
locity oscillations of stars with exoplanets, at a resolv- 
ing power of 42 000 in the wavelength range 3850-6800 A. 
The stability of the system makes it possible to compare 
spectra taken at different epochs. Spectrum extraction, 
wavelength calibration and radial velocity measurement 
by cross-correlation have been performed at the telescope 
with the on-line data reduction software. 

The current version of the Elodie library includes 1962 
spectra available in the Hypercat database 1 , most of the 
spectra having a signal to noise ration (S/N) at 550 nm 
greater than 100. We have selected 208 spectra of G 
dwarfs with the following criteria : 0.55 < B — V < 0.75 
{{B - V) B ~ 0.65) and 4 < M v < 5.6 (M V0 ~ 4.82). 
Absolute magnitudes have been computed from Hipparcos 
parallaxes, the selected stars having relative errors smaller 
than 10%. The list includes 8 spectra of the Sun (Table 

m 



Table 1. List of the solar spectra used in this study. 



Hypercat 


date of 


object 


FWHM 


S/N at 


number 


observation 




km.s -1 


550 nm 


00903 


14/01/1998 


Moon 


11.061 


381.4 


00904 


22/12/1999 


Moon 


11.050 


268.5 


00905 


22/12/1999 


Moon 


11.050 


139.6 


00906 


22/12/1999 


Moon 


11.061 


156.5 


00907 


27/03/2000 


Ceres 


11.017 


117.5 


00908 


24/01/2000 


Moon 


11.057 


200.0 


00909 


24/01/2000 


Moon 


11.054 


224.9 


01964 


22/08/2000 


Moon 


12.126 


404.3 



1 www-obs.univ-lyonl.fr/hypercat /ll/ spectrophotometry. 



The stellar spectra were compared to solar ones 
with the TGMET code developed by Katz et al fiMSjl . 
TGMET is a minimum distance method which measures 
similarities between spectra in a quantitative way. The 
TGMET comparison between 2 spectra includes the fol- 
lowing steps : 

— straightening of each order 

— removal of bad pixels, cosmic hits and telluric lines 

— wavelength adjustment 

— mean flux adjustment by weighted least squares, order 
by order 

The wavelength adjustment shifts the comparison 
spectrum to the radial velocity of the solar spectrum and 
resamples it to the same wavelenghts. It implies an inter- 
polation between wavelengths which is performed with the 
quadratic Bessel formula. The flux fitting of the compar- 
ison spectrum to the solar spectrum assumes that the 2 
spectra differ by a simple factor (the 2 stars having roughly 
the same temperature it is not necessary to introduce a 
slope) . Once the two spectra have been put on a common 
scale, a distance between them can then be computed. As 
explained in Katz et al (1998), instead of adopting the 
real reduced \ 2 of the fit as the distance between 2 spec- 
tra, which would imply taking into account the noise on 
each pixel, the response curve was chosen as the weighting 
function. This distance was adopted after many tests and 
was proven to produce the most satisfactory results, espe- 
cially at high S/N. Its advantage is that it gives a similar 
weight to the continuum and to the wings and bottom of 
absorption lines, contrary to a weighting function based 
on the photon noise. 

Katz et al's algorithm includes a convolution step 
which is not included in the present work. A convolution 
should be performed in order to put the 2 compared spec- 
tra at exactly the same resolution. A difference in resolu- 
tion between two spectra is the result either of a variation 
of the instrumental resolution between the two exposures 
or of intrinsic physical properties of the observed stars 
like rotation, macroturbulence or binarity which enlarge 
their spectral profile. But as we are looking for solar twins, 
these intrinsic properties are important in the criterion of 
similarity and should not be erased. Moreover ELODIE is 
a very stable instrument and its resolution does not vary 
significantly with time. It can be seen in TableQ]that all of 
our solar spectra have a resolution of 11 km.s -1 (FWHM), 
except the spectrum 01964 which is sligthly degraded (12 
km.s -1 ). These considerations led us to supress the con- 
volution step in the TGMET algorithm. 

In practice we have limited the comparison to the 
wavelength range 4400-6800 A (orders 21 to 67) and elim- 
inated the under- illuminated edges of the orders. Finally 
distances between spectra have been computed over nearly 
32000 wavelengths. Table gives an example of the 
TGMET output, for solar spectra 00903 (S/N=381.4) and 
00907 (S/N=117.5). The output in the two cases is consis- 
tent, with however some differences : HD088072 is within 
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the 20 closest neigbours of 00903 but not of 00907, the 
opposite is the case for HD071148 and HD042618. These 
differences, probably related to observing conditions, are 
smoothed when combining the TGMET results obtained 
for the 8 solar spectra, the combination being performed 
by averaging the distances, order by order, giving a dif- 
ferent weight to several orders (see below). The score ob- 
tained by each solar analog was our criterion of closeness 
to the Sun, leading the final Top Ten list : HD146233, 
HD168009, HD010307, HD089269, HD047309, HD095128, 
HD042618, HD071148, HD186427, HD010145. 

Fig. U] shows the fit of HD146233 to spectrum 00903, 
for order 39 including the Mgl triplet and order 64 in- 
cluding the H a line. For order 39, the fit has been per- 
formed on 773 points; the mean difference between the 
solar flux and the fitted fluxes of HD146233 is 3 electrons 
with a standard deviation of 1157 electrons, corresponding 
to 1% of the mean flux (108772 electrons). For order 64, 
the difference is also 1% (743 points, mean difference of 
- 57 electrons, standard deviation of 2166 electrons, mean 
flux of 210319 electrons). 

Fig. shows for the Top Ten solar analogs their dis- 
tance to the solar spectrum 00903, order by order. It is 
very clear from Fig. [21 that HD146233 is very similar to 
the Sun and that HD168009 is not very far behind. The 
closeness of these 2 stars is confirmed for the 7 other so- 
lar spectra. The largest discrepancies occur for order 63 
(648- 652.5 nm). The examination of this order indicates 
that it is polluted by telluric lines which were not com- 
pletely removed. Telluric lines unfortunatly affect also or- 
der 64 which is our best indicator of temperature thanks 
to the H a line. The dispersion obtained on order 64 is 
much higher than that on order 31 which includes the Hp 
line. However, the Hp line being at the edge of the order it 
has a lower weight in the fit because of undcr-illumination. 
A large dispersion on order 39, which includes strong fea- 
tures due to the Mgl triplet, is also seen for each solar 
spectrum. This region is known to be very sensitive to the 
3 atmospheric parameters (T e g, logg, [Fe/H]) and con- 
sequently powerful for discriminating solar twins. These 
considerations led us to adopt a higher weight of 3 on or- 
der 39, a half weight on order 64 and a null weight on order 
63 when combining the information on all the orders. 

3. Atmospheric parameters and photometry 

We report in this section colour indexes available in the 
literature for our Top Ten solar analogs and visual ab- 
solute magnitudes deduced from their Hipparcos paral- 
laxes (Table 0}. We also review recent determinations of 
their atmospheric parameters. These data are compared to 
those of the Sun, and uncertainties which affect them are 
discussed. Finally we report their occurrence in previous 
studies of solar analogs. 

The B — V and U — B colours come from the General 
Catalogue of Photometric Data (Mermilliod et al 11997(1 . 
except for HD047309 for which we have taken the B — V 
colour from Tycho2 (Hog et al 2000} transformed to the 



Johnson system. The b-y colours are extracted from the 
catalogue by Hauck & Mermilliod (|1998[l . 

The photometry of the Sun is a matter of debate. 
Neckel i|T9l^ has determined (B - V) & = 0.650 ± 
0.005, (U - B) e = 0.195 ± 0.005, M v © = 4.82 ± 0.025, 
values which are adopted as basic solar data in Allen's 
Astrophysical Quantities (2000). Cayrel de Strobel (JT996 ) 
gives a compilation of solar (B — V) colours measured 
by different techniques and determines from the relations 
colour vs Teff : (B - V) e = 0.642 ± 0.004 , (6 - y) e = 
0.404±0.005. Porto de Mello & da Silva l(TM7| obtain with 
a similar method : (B-V) e = 0.648 ± 0.006^ (U - B) Q = 
0.1 78 ±0.013. We list in Table H a reasonable range of 
values for the Sun's colours and absolute magnitudes. 



3.1. HD 146233 

HD146233 (18 Sco) was adopted at the Solar Analogs 
workshop at Lowell Observatory as one of the best so- 
lar twins. Our study confirms with independent data and 
methods the result of Porto de Mello & da Silva l)1997[> 
quoting HD146233 as THE closest ever solar twin. In the 
optical range, its spectrum is indistinguishable from that 
of the Sun (Fig. P). Before that, Hardorp l|1978f) using 
UV spectrophotometry mentioned this star as a solar ana- 
log but with the comment "spectrum similar to solar, 
some absorption features weaker". However, this study 
was based on a single low-resolution spectrum. This dis- 
crepancy is discussed by Porto de Mello & da Silva Q1997[) 
and by Cayrel de Strobel ifEMj) . 

Only 2 determinations of atmospheric parameters are 
available for this star, in very good agreement, giving so- 
lar values within the error bars. HD 146233 seems to be 
more luminous than the Sun by 0.05 mag. Its parallax is 
very accurate (ir — 71.30 ± 0.89 mas), but one may le- 
gitimely wonder if its photometry is sufficiently accurate 
to consider this excess of luminosity real. As a matter 
of fact HD146233 is part of the Catalogue of suspected 
variables by Kukarkin et al Ijl981|) who found a possible 
amplitude of 0.11 mag. V magnitudes, measured by sev- 
eral authors between 1957 and 1978 and available in the 
GCPD (Mermilliod et al 1997), range effectively between 
V=5.48 and V=5.56. The average V = 5.504 ± 0.015 was 
used to compute an absolute magnitude of My — 4.77. 
Thus if HD146233 is slightly variable, a higher luminosity 
than the Sun cannot be clearly established. But more re- 
cently HD146233 was identified by Adelman pToTj) to be 
part of the 681 most photometrically stable stars during 
the 5 years of the Hipparcos mission, with an amplitude of 
0.01 mag. Its median magnitude in the Hipparcos system 
is Hp — 5.6265 ± 0.0005. According to the photometric 
transformation calibrated by Harmanec (|1998f) . the corre- 
sponding apparent visual magnitude is V — 5.493 ± 0.003 
which confirms its higher luminosity. 
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Fig. 1. Comparison of one of the solar spectra (dots) with HD146233 (continuous line) 
Mgl triplet and H a line. 
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Table 2. The 20 closest spectra of the solar spectra 00903 and 00907, deduced by the TGMET code. 



n° 


star 


S/N 


distance 


n° 


star 


S/N 


distance 




00903 






00907 




01964 


Sun 


404.3 


2.220 


00903 


Sun 


381.4 


1.158 


00909 


Sun 


224.9 


2.332 


00909 


Sun 


224.9 


1.175 


00908 


Sun 


200.0 


2.410 


01964 


Sun 


404.3 


1.183 


00490 


HD146233 


236.8 


2.789 


00908 


Sun 


200.0 


1.188 


00906 


Sun 


156.5 


2.910 


00490 


HD 146233 


236.8 


1.243 


00905 


Sun 


139.6 


3.110 


00906 


Sun 


156.5 


1.250 


01633 


HD168009 


204.9 


3.305 


00905 


Sun 


139.6 


1.275 


00904 


Sun 


268.5 


3.392 


01633 


HD168009 


204.9 


1.329 


00907 


Sun 


117.5 


3.474 


00904 


Sun 


268.5 


1.382 


01187 


HD047309 


119.7 


3.841 


00039 


HD010307 


198.4 


1.457 


01188 


HD047309 


108.8 


3.864 


01634 


HD168009 


134.6 


1.458 


01634 


HD168009 


134.6 


3.875 


01187 


HD047309 


119.7 


1.471 


00039 


HD010307 


198.4 


3.937 


00346 


HD071148 


117.0 


1.472 


00895 


HD089269 


225.9 


3.978 


00895 


HD089269 


225.9 


1.478 


00258 


HD047309 


100.9 


4.176 


01188 


HD047309 


108.8 


1.493 


00387 


HD088072 


86.5 


4.236 


00400 


HD095128 


181.9 


1.555 


00699 


HD186427 


139.9 


4.283 


00258 


HD047309 


100.9 


1.557 


00400 


HD095128 


181.9 


4.336 


00981 


HD010307 


162.9 


1.560 


00981 


HD010307 


162.9 


4.342 


01125 


HD042618 


132.5 


1.591 


00038 


HD010145 


153.2 


4.346 


00699 


HD186427 


139.9 


1.594 



3.2. HD 168009 

HD 168009 has been quite well studied but has never been 
mentioned as a good solar analog, despite its being part 
of the list of " bright stars with UBV-colours close to those 
of the Sun" established by Neckcl (1986). According to its 



spectroscopic gravity and absolute magnitude, HD168009 
is more luminous than the Sun. Its absolute magnitude 
My = 4.52 is quite reliable and relies on a parallax of 
7r = 44.08 ± 0.51 mas. Several estimations of its appar- 
ent visual magnitude are in good agreement : V=6.295 
according to the GCPD, V=6.309 according to Simbad, 




order 



+ HD1 46233 *HD1 68009 °HD0 1 0307 x HD089269 D HD047309 
A HD0951 28®HD04261 8°HD071 1 48^HD186427^HD010145 

Fig. 2. Distance of the Top Ten solar analogs to solar spectrum 00903, order by order. 



V=6.307 according to Tycho2. Its B — V colour index is 
slightly more uncertain : B — V = 0.635 according to the 
GCPD, B-V = 0.604 according to Simbad, B—V = 0.646 
according to Tycho2, but suggests however a higher tem- 
perature than the Sun. This is confirmed by several recent 
estimates of Teff available in the literature and spanning 
values from 5719K (Chen et al l2UUU|) to 5833K (Blackwell 
& Lvnas-Grav I1998|) with a mean value of 5801K. This 
large dispersion is a good illustration of the lack of a 
common temperature scale, even for bright nearby stars. 
HD 168009 is also part of a catalogue of high precision near 
infrared photometry by Kidger & Martin-Luis l|2003|l who 
give J=5.133, H=4.840, K=4.783, values which differ by 
less than 0.005 mag from those measured by Alonso et al 
ifEMjl . A colour index V-K=1.512 leads to Teff=57 30 K 
with the relation established by Alonso et al (1996b). 

The abundance of several elements (O, Na, Mg, Al, Si, 
K, Ca, Ti, V, Cr, Ni, Ba) have been measured by Chen 
et al (2000) to be solar within the error bars. Ba, Eu and 
Sr abundances have also been measured by Mashonkina 
& Gehren (2001), leading to the same conclusion. 



et al l|1992[) . the low mass companion being 1000 times 
fainter. According to a detailed analysis of Hipparcos data 
by Martin et al (19!JH|) the system has a total mass of 
(0.931 ± 0.178)M Q , a primary mass of (0.795 ± O.159)M 
and a secondary mass of (0.136±0.053)M Q . At a distance 
of 12.6 pc (12.4 pc when the binarity is considered), it is 
the nearest of our Top Ten solar analogs. It is a well stud- 
ied star with many measurements in good agreement of 
its apparent visual magnitude and B ~ V colour, and it is 
part of the catalogue of the least variable stars compiled 
by Adelman l|2001fl despite its binarity. We notice that 
the temperature given by Chen et al is similar to that of 
the Sun but significantly lower than given by other au- 
thors, as was also the case for HD168009. Hardorp (1978) 
mentioned HD010307 with the comment "some aborption 
features appreciably weaker than solar" which is in agree- 
ment with a higher temperature. 

Cayrel de Strobel l|1996j> and Fesenko l|1994[> have in- 
cluded HD010307 in their list of solar analogs but not with 
a high rank. It was only mentioned to "deserve study" in 
the conclusions of the Lowell Workshop. 



3.3. HD010307 

Like HD168009, HD010307 seems to be hotter and more 
luminous than the Sun. Allende Prieto et al l|1999|l quote a 
mass of 0.94Mq and an astrometric gravity of logg=4.29, 
in very good agreement with the averaged spectroscopic 
gravity of logg=4.26 quoted in Table H HD010307 is 
in fact a binary system which was resolved by Henry 



3.4. HD089269 

Very few papers mentioning HD089269 are available in 
the literature. It was never recorded as a solar analog de- 
spite its colour index B — V = 0.654 being similar to that 
of the Sun. Its visual magnitude V=6.633 combined to 
its trigonometric parallax ir = 48.45 ± 0.85 mas leads to 
an absolute magnitude My = 5.06 indicating that it is 
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less massive than the Sun. Mishcnina et al (2003) find it 
100K colder than the Sun and significantly more metal 
poor. The good score obtained with TGMET indicates 
that, globally, the combined effects of temperature and 
iron abundance may give similar absorption features as in 
the Sun. 

3.5. HD047309 

HD047309 is unknown as a solar analog. The only data 
are Stromgren photometry, and data from Hipparcos and 
Tycho2. Contrary to HD089269, it seems to be slightly 
hotter, more luminous and metal-rich. It is also the most 
distant of our sample at 42.4 pc. 

3.6. HD095128 

We come back to well a known star with HD095128 (47 
UMa) known to have two giant planets orbiting around it. 
Consequently it has been very well studied. Its tempera- 
ture is significantly higher than that of the Sun. Again 
the temperature given by Chen et al is the lowest of 
the list with 5731K whereas Santos et al (|2003p find the 
highest temperature, with 5925K. The 2 colour indexes 
B — V = 0.606 and b — y = 0.391 confirm a higher temper- 
ature than the Sun. Despite the large dispersion in Teff, 
the dispersion in [Fe/H] is low, with an average exactly 
solar. HD095128 is more luminous than the Sun. Allende 
Prieto et al (1993) quote a mass of 0.96M Q and an as- 
trometric gravity of logg=4.23, in good agreement with 
the averaged spectroscopic gravity of logg=4.28 quoted in 
Table^J It shows that HD095128 has already evolved from 
the main sequence. 

3.7. HD042618 

The only mention of HD042618 is by Fesenko (flUM)l in 
his solar type star study. The 2 available estimates of its 
temperature differ by 120K, but B — V = 0.632 suggests 
that the hottest one, 5775K, is the most likely. The cold 
temperature scale adopted by Reddy et al ( 2002 ) implies a 
low metallicity of [Fe/H]=-0.16, which might be closer to 
the solar one in fact. It is not very clear however, because 
its absolute magnitude My = 5.05 combined with a solar 
temperature would not be compatible with being on the 
solar composition ZAMS. 

3.8. HD071148 

Like HD042618, HD071148 was recorded by Fesenko 
(1994) but did not receive much attention as a solar ana- 
log. Its effective temperature is subject to a controversy 
between partisans of a low temperature scale (Reddy et al 
2002 and Chen et al l2003|) giving a temperature of about 
5710K and Mishenina et al §003) giving a temperature 
140K higher. Again we find that the colour indexes in- 
dicate that this star is probably slightly hotter than the 



Sun. The 3 authors agree on the fact that its metallicity 
is nearly solar. 

Interestingly, the radial velocity of HD071148 has been 
monitored during several years by Naef et al ( 2003 ) who 
found that it is constant within 10 ins -1 , ruling out the 
presence of a low mass companion. 

3.9. HD 186427 

HD186427 (16 Cyg B) is one of the best solar twin candi- 
dates of the Lowell Fall Workshop and also a planet-host 
star. Its UV spectrum was qualified as "indistinguishable 
from solar" by Hardorp (1978), a similarity confirmed by 
Fernley et al (1996|. There is a remarkable agreement be- 
tween the 10 authors who have estimated its temperature, 
colder than the Sun by only ^20K. Only Laws & Gonzalez 
( 2001 ) and Gonzalez ( 1998 ) used a significantly lower tem- 
perature scale. HD186427 differs from the Sun mainly by 
its higher metal content, and by a higher luminosity. 

3.10. HD010145 

HD010145 has never been mentioned as a solar analog and 
was little studied before. The recent spectroscopic analysis 
performed by Mishenina et al ( 2003 ) shows that it is colder 
than the Sun by 100K, but with a similar gravity and 
metal content. It is the coldest of our Top Ten together 
with HD089269, the latter one having bluer colour index 
consistent with its lower metallicity. HD010145 has a lower 
luminosity than the Sun, indicating either a lower age or 
a lower mass. 

4. Li content, activity, ages 

In this section we compare qualitatively the Li content of 
our Top Ten solar analogs with that of the Sun. The solar 
photosphere is known to be highly depleted in Li, as is 
the case for many solar type stars. This depletion is how- 
ever subject to various interpretations, involving rotation, 
convection, or the presence of a planetary system. The 
correlation of age with Li depletion has also been often 
discussed but not fully established. 

The 7 Li doublet resonance lines at 670.78 nm and 
670.79 nm are well placed in the middle of the 66th order 
of the ELODIE spectra. In Fig.0 for each of our Top Ten, 
the Li region was superposed on the solar spectrum, show- 
ing that 6 of them are depleted in Li like the Sun whereas 
the 4 others show a pronounced Li feature. The most pro- 
nounced feature concerns HD071148 and HD010307. A 
weaker feature is also present in the spectra of HD095128 
and HD146233. 

More clearly than the Li content, the chromospheric 
activity of a solar type star is directly connected to its age. 
Thus it would have been extremely interesting to look at 
the chromospheric activity revealed by the central depth 
of the Ca II H and K lines, and of the Ca II triplet lines 
at 852 nm. Unfortunatly the NIR Ca II triplet is not in 
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the spectral range of ELODIE and the H and K lines ap- 
pear on the border of the 2nd and 3rd orders which are 
underilluminated. The core of the H a line is also an indi- 
cator of chromospheric activity, but we were not able to 
detect any significant difference of depth, even for the 4 
stars having a higher Li content. According to Soderblom 
(1985). HD071148, HD010307 and HD095128 show Call 
emission strengths and rotation similar to the Sun, sug- 
gesting a weak activity. Hall & Lockwood l|2000|) found 
an activity cycle in HD146233 with an amplitude slightly 
greater than that of the Sun. Thus these 4 stars present- 
ing a pronounced Li feature do not seem to be enormously 
more active than the Sun. Two other analogs are part of 
Soderblom's study, HD010145 and HD168009; they do not 
exhibit evidence of a higher activity than the Sun. 

Several of our Top Ten had their age estimated by 
Ibukiyama & Arimoto <|2TjH2)> : 7.02 Gyr for HD168009, 
7.32 Gyr for HD010307, 6.92 Gyr for HD095128, 6.65 Gyr 
for HD071148. According to Cayrel de Strobel & Friel 
(JTnUBI), HD146233 has an age of 6 Gyr. Thus one can find 
stars older than the Sun which have a higher Li content. 

5. Solar analogs selected by other methods 

The ELODIE library includes the spectra of several stars 
which were considered in previous studies as high rank 
solar analogs. We will focus here on HD 2 170 14 (51 Peg), 
HD028099 (Hy VB 64) and HD186408 (16 Cyg A). 

With our purely differential method HD 186408 is 
within the 20 best solar analogs among our list of 208 
dwarfs. A dozen recent determinations of its atmospheric 
parameters are available and give on average Teff=5780K, 
logg=4.26, [Fe/H]=+0.07. Its visual absolute magnitude is 
My = 4.29, to be compared to My = 4.56 for its compan- 
ion HD186427 and M VQ = 4.82 showing that HD186408 
is more massive and evolved. 

HD217014 and HD028099 obtained a lower score in 
the TGMET output despite their similar effective tem- 
perature and gravity to the Sun. These 2 stars have also 
been very much studied, especially HD217014 because of 
its planet. The latest determination of its atmospheric 
parameters le d to T eff=5805K, logg=4.51, [Fe/H]=+0.21 
(Santos et al 2003). Thus stronger metallic lines explain 
that it is not a good spectroscopic analog of the Sun. 
HD028099 was found to have Teff=5800K, logg=4.40, 
[Fe/H]=+0.10 by Paulson et al <j27)TI3|l . This star is not 
considerably more metal rich than the Sun, but it is known 
to be younger and to have a high chromospheric activity. 
Looking closely at its spectrum shows that its lines are 
not as deep as those in the solar spectrum. 

We have also searched for photometric solar analogs 
in the ELODIE library by selecting stars having the same 
colours and absolute magnitude as the Sun. Four stars 
fall in the narrow range 4.6 < My < 5.0, 0.63 < B - V < 
0.68,0.16 < U-B < 0.23,0.40 < b-y < 0.42 : HD001835 
(BE Cet), HD076151, HD146233, HD159222. HD076151 
in not in our Top Ten but it is at the 14th position. It 
is a well studied star, with recent determinations of its 
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HD0 


89269 



670.4 670. 




Fig. 3. Comparison between the Li I 670.7 nm region of 
the Sun (spectrum 00903, dotted line) and of each Top 
Ten solar analog 



atmospheric parameters in good agreement, giving on av- 
erage Teff=5774K, logg=4.39, [Fe/H]=+0.06. It was also 
one of the good solar twin candidates discussed by Cayrel 
de Strobel : 199(i;. However, this star is younger than the 
Sun (3.04 Gyr estimated by Ibukiyama & Arimoto 2002), 
with a stronger activity and faster rotation (Pizzolato et 
al 2003 ) and presents a pronounced Li feature. In combi- 
nation with a higher metal abundance, this may explain 
why it was not at higher rank in the TGMET output. 
HD001835 is a variable star of the BY Dra type which 
renders its photometry suspect. It was also discussed by 
Cayrel de Strobel (1996) as a good solar analog for the 
temperature and gravity but not for Li, chromospheric 
activity and age. Its young age is assessed by its member- 
ship of the Hyades moving group and its high activity is 
confirmed by Pizzolato et al ((2)03). It is also more metal- 
rich than the Sun ([Fe/H]=+0.13 estimated by Mishenina 
et al 120031) and its distance to the solar spectra computed 
by TGMET is very large. Its activity is clearly seen in its 
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H a line, which is shallower than in the Sun. HD159222 
the 11th star in the TGMET output. It is thus a good 
photometric and spectroscopic solar analog. Moreover its 
age is also very similar (4.56 Gyr estimated by Ibukiyama 
& Arimoto 2002 ) . Several determinations of Teff are avail- 
able for this star, showing an impressive dispersion : 5708 
K by Blackw ell & Lynas-Gray lfTM5|l . 577 K b y Alonso 
et al <ll996at . 58 34 K b y Mishenina et al lIMiSTl . 5852 K 
by di Benedetto fTOSBI . 

We have also observed that many spectra were pol- 
luted by telluric lines in the red orders, and that order 
39 including the Mgl triplet was a powerful discrimina- 
tor of solar resemblance. We have thus performed the 
TGMET comparison of the 8 solar spectra with the li- 
brary using only this order. This has greatly modified 
the order of our list, HD146233 keeping however its high- 
est rank. The five closest stars are part of the Top Ten 
list : HD146233, HD047309, HD168009, HD042618 and 
HD186427. The 4 stars of the Top Ten list with the lowest 
temperatures (HD010145 and HD089269) and the high- 
est temperatures (HD010307 and HD095128) are pushed 
away. Five new stars appear in the 10 closest solar analogs 
: HD195034, HD159222 (also a good photometric analog), 
HD187123, HD186104, HD005294. The photometric ana- 
log HD076151 is at the 11th position. 

The parameters of the solar analogs discussed in 
this section are listed in Table |21 All of them, except 
HD005294, are more metal-rich than the Sun. 

6. Discussion 

It is interesting to note that, despite the great similitude 
of the optical spectrum of our Top Ten solar analogs to 
that of the Sun, their atmospheric parameters can dif- 
fer significantly. Effective temperatures span ±100K on 
both sides of the solar value, logg values span the interval 
[4.09 ; 4.58] and [Fe/H] span the interval [-0.23 ; +0.11]. 
Several interpretations can explain this dispersion. On the 
one hand, authors do not use the same temperature scale 
and model atmospheres. Temperature scales can differ by 
more than 150K as we have seen. It is very important 
that authors agree on temperatures because this param- 
eter has a strong impact on abundance determinations. 
Uncertainties which affect the stellar parameter determi- 
nations have been discussed by Asplund (J2D03 ) to be of 
the order of 50K to 100K in Teff, 0.2 dex for logg and 0.1 
dex for [Fe/H]. On the other hand, we cannot expect find- 
ing a perfect twin having all its parameters exactly solar, 
especially in an incomplete sample. Finally, temperature 
and metallicity have contrary effects on the overall spec- 
trum which may compensate in some cases (ex HD089269 
or HD187123). It is also possible that other effects act on 
the spectra. Observing conditions and telluric lines are the 
most obvious, but intrinsic stellar properties also have an 
influence on the spectrum. We have mentioned chromo- 
spheric activity and rotation, but the abundance of other 
elements than iron, turbulent motions, spots on the stel- 
lar surface may be different than in the Sun. For instance, 



when using TGMET only in the Mgl triplet region, the 
Mg abundance of the star may have a strong weight. 

Colour indexes of the Top Ten solar analogs also span 
intervals as large as 0.073 in B — V, 0.032 in 6 — y, 0.099 in 
U — B, and absolute magnitudes range from M v = 4.31 to 
M v = 5.06. Fig. ^represents their distribution in Teff vs 
colour and HR diagrams together with the other analogs 
discussed in Sect.[S]and the rest of the ELODIE library, re- 
stricted to -0.25 < [Fe/H] < +0.15. The Stromgren index 
b — y is clearly the one which presents the lowest disper- 
sion in its relation with Teff. Like atmospheric parameters, 
magnitudes and colours are affected by uncertainties and 
a lack of homogeneity. Absolute magnitudes are computed 
from excellent parallaxes but averages of old and inhomo- 
geneous apparent magnitudes. Tycho2 (H0g et al l2000J) is 
a recent photometric catalogue of good quality but its B 
and V passbands do not correspond to the Johnson stan- 
dard system and transformations, also affected by calibra- 
tion uncertainties, have to be used. A small fraction of the 
dispersion may also be due to interstellar absorption, even 
if our targets are closer than 50 pc. But we interpret the 
larger part of the dispersion to mean that our incomplete 
sample of 208 G dwarfs includes stars of various ages and 
states of evolution resulting in a variety of astrophysical 
properties, because the ELODIE library was built to rep- 
resent the stellar content in the solar neighbourhood, not 
solar analogs. 

We have seen in previous sections that good photo- 
metric analogs are not always good spectroscopic analogs. 
HD001835 is a good example of a star having similar 
colours and absolute magnitude as the Sun but which is 
considerably different in respect to its age, activity and 
metal content. Thus photometry is not able to discrimi- 
nate between these effects whereas high resolution spec- 
troscopy can. In contrast the direct comparison of high 
resolution spectra used alone classify as solar analogs stars 
having a large range of atmospheric parameters. We con- 
clude that a good strategy to find other solar twins than 
HD146233, and perhaps better ones, would be to select 
photometric analogs in large catalogues, then select with 
Hipparcos those that have a similar absolute magnitude 
to the Sun, then submit their high resolution spectrum to 
the TGMET comparison. We have scanned the GCPD, 
Stromgren and Hipparcos catalogues with the drastic cri- 
terion used in Sect. [5] to identify photometric analogs and 
found only 27 candidates, 4 being already in the ELODIE 
library, 15 others observable with ELODIE. We plan to 
observe them soon in order to complete this work. 

7. Conclusion 

We have presented the 10 stars of the ELODIE library 
which exhibit the closest optical spectrum to the Sun at 
a resolution of 42000. They have colours, absolute mag- 
nitudes, atmospheric parameters and Li content which 
span a range of values larger than expected. It is surpris- 
ing that a star like HD089269, colder, more metal poor 
and less luminous than the Sun is at the 4th position, 
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Table 3. Fondamental data for solar analogs mentioned in Sect|SJ Only the most recent determination of atmospheric 
parameters is listed. The second colomn indicates how the star was selected as a solar analog (P : photometry, L : 
Lowell workshop, Mg : TGMET on the Mg I triplet region) 



Star name 


method 


B - V 


b-y 


U - B 


M v 


Teff 


log g 


[Fe/H] 


source 


HD186408 


L 


0.645 


0.410 


0.187 


4.29 


5803 


4.20 


+0.02 


Mishenina et al ( 2003 ) 


HD217014 


L 


0.665 


0.416 


0.224 


4.53 


5805 


4.51 


+0.21 


Santos et al <I2003|I 


HD028099 


L 


0.660 


0.411 




4.75 


5800 


4.40 


+0.10 


Paulson et al 


HD076151 


P 


0.662 


0.413 


0.217 


4.83 


5776 


4.40 


+0.05 


Mishenina et al H2003I 


HD001835 


P 


0.659 


0.409 


0.226 


4.84 


5790 


4.50 


+0.13 


Mishenina et al (I2003H 


HD005294 


Mg 


0.650 


0.401 


0.174 


5.03 


5779 


4.10 


-0.17 


Mishenina et al <|2003fl 


HD159222 


P+Mg 


0.637 


0.406 


0.172 


4.67 


5834 


4.30 


+0.06 


Mishenina et al (I2003H 


HD186104 


Mg 


0.631 


0.412 




4.62 


5753 


4.20 


+0.05 


Mishenina et al (1200311 


HD187123 


Mg 


0.619 


0.405 




4.43 


5855 


4.48 


+0.14 


Santos et al ||2QQ3j 


HD195034 


Mg 


0.610 


0.408 




4.84 

















? 






0.6 


0.55 

a-v 


0.7 




Fig. 4. Distribution of the ELODIE library (small dots), 
the Top Ten solar analogs (stars), other solar analogs 
(Sect. diamonds) and the Sun in Teff vs colour and 
HR diagrams. The solar twin HD146233 is shown with a 
larger symbol. 



whereas HD076151, having similar colours, absolute mag- 
nitude and atmospheric parameters is only at the 14th 
position. Activity may play an important role in discrim- 
ination. Wc have shown for instance that the good pho- 
tometric analog HD001835 was a very bad spectroscopic 
analog because of its high activity. One also has to take 
into account, when comparing colours, absolute magni- 
tudes and atmospheric parameters to those of the Sun, 
that these quantities are affected by significant uncertain- 
ties. Effective temperatures are particularly in question, 
with determinations for the same star differing by nearly 
200K in some cases. Our method consisting in measur- 
ing distances between spectra is powerful but it is also 
affected by uncertainties due to observing conditions, es- 
pecially the pollution by telluric lines, which may perturb 
the order of the classification. 

Among our Top Ten, several stars have never been 
mentioned before as solar analogs and have been very lit- 
tle studied. They are good candidates for planet hunt- 



ing, especially HD047309 which is slightly more metal rich 
than the Sun. Two of our solar analogs, HD095128 and 
HD186427, are already known to have planets. HD159222 
and HD076151 are also good candidates because they are 
good spectroscopic analogs (in the Top 15) and good pho- 
tometric analogs. 

The conclusion of this work is that none of the methods 
to search for solar twins is satisfactory when used by itself. 
The methods that have been already used are the com- 
parison of colour indexes, of absolute magnitudes, of UV 
spectral energy distributions, of atmospheric parameters 
and of high resolution optical spectra. All these methods 
are affected by uncertainties and none of them is able to 
describe sufficiently all the stellar properties. Combining 
them is the best way to minimize their drawbacks, uncer- 
tainties and insufficiencies. Finally HD146233 is the only 
star in the ELODIE library which merits the title of solar 
twin because it has passed the filter of all methods. It is 
not however a perfect twin and differs from the Sun by 
its higher Li content, slightly higher age (6 Gyr against 
4.6 Gyr for the Sun) and higher luminosity (My = 4.77 
against M V0 = 4.82). 
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Table 4. Fondamental data for our Top Ten solar analogs. Our adopted atmospheric parameters in bold characters are the mean values from the literature. 
The second column is the averaged TGMET distance of the corresponding star to the 8 solar spectra. 



Star name 


TGMET 


B — V 


b-y 


U-B 


distance 


Absolute 


Teff 


log g 


[Fc/H] 


source 




score 








in pc 


magnitude 








(Teff, logg, [Fe/H]) 


Sun 




0.64-0.66 


0.40-0.41 


0.17-0.20 





4.80-4.84 


5777 


4.44 


0.00 


- 


HD 146233 


2.019 


0.651 


0.401 


0.174 


14.0 


4.77 


5799 


4.40 


+0.01 


Mishenina et al P003 1 
















5789 


4.49 


+0.05 


Porto de Mello & da Silva fTWl 
















5794 


4.44 


+0.03 




HD168009 


2.195 


0.635 


0.410 


0.115 


22.7 


4.52 


5826 


4.10 


-0.01 


Mishenina et al ( 2003 1 
















5719 


4.08 


-0.07 


Chen et al P5oTJ|l 
















5826 


- 


- 


di Benedetto 
















5833 


- 


- 


Blackwell & Lynas-Gray fT558l 
















5781 


- 


- 


Alonso et al I1996al 
















5801 


4.09 


-0.04 




HD010307 


2.516 


0.616 


0.389 


0.113 


12.6 


4.45 


5881 


4.30 


+0.02 


Mishenina et al ( 2003 1 
















5776 


4.13 


-0.05 


Chen et al J5oTO 
















5825 


4.33 


-0.04 


Gratton et al 1199611 
















5874 


- 


- 


Alonso et al l|1996a|l 
















5898 


4.31 


-0.02 


Edvardsson et al 115531 
















5848 


4.26 


-0.02 




HD089269 


2.602 


0.654 


0.420 


0.156 


20.6 


5.06 


5674 


4.40 


-0.23 


Mishenina et al p003 1 
















5674 


4.40 


-0.23 




HD047309 


2.628 


0.623 


0.412 


- 


42.4 


4.47 


5791 


- 


- 


Kovtyukh et al ( 2003 ) 
















5791 


4.40 


+0.11 


Mishenina priv. com. 
















5791 


4.40 


+0.11 




HD095128 


2.758 


0.606 


0.391 


0.126 


14.1 


4.31 


5887 


4.30 


+0.01 


Mishenina et al ( 2003 1 
















5861 


4.29 


+0.05 


Laws et al ( 20031 
















5925 


4.45 


+0.05 


Santos et al (120031 
















5788 


4.31 


+0.03 


Zhao et al (121)521 
















5731 


4.16 


-0.12 


Chen et al (1211501 
















5892 


4.27 


0.00 


Zhao & Gehren (1277501 
















5892 


4.27 


0.00 


Fuhrmann et al (119981 
















5800 


4.25 


+0.01 


Gonzalez (1T5551 
















5892 


4.27 


0.00 


Fuhrmann et al (119971 
















5811 


4.09 


0.00 


Gratton et al Ijl996|l 
















5882 


4.34 


+0.01 


Edvardsson et al 1)1993^ 
















5860 


4.28 


0.00 





CO 

o 

g" 

e 
p 

B 
e 



CO 

o 



P 

2L 

o 



Star name 


TGMET B — V 


b-y 


U-B 


distance 


Absolute 


Tcff 


log g 


[Fe/H] 


source 




score 






in pc 


magnitude 








(Teff, logg, [Fe/H]) 


Sun 


0.64-0.66 


0.40-0.41 


0.17-0.20 





4.80-4.84 


5777 


4.44 


0.00 


- 


HD042618 


2.857 0.632 


0.404 


0.127 


23.1 


5.05 


5775 


- 


- 


Kovtyukh ct al ( 2003 ) 














5653 


4.58 


-0.16 


Reddy et al ll2"UU2l 














5714 


4.58 


-0.16 




HD071148 


2.907 0.625 


0.399 


0.120 


21.8 


4.65 


5850 


4.25 


0.00 


Mishenina et al ( 2003 ) 














5716 


4.34 


+0.03 


Chen et al (I3U031 














5703 


4.46 


-0.08 


Reddy et al ll2"UU2l 














5756 


4.35 


-0.02 




HD186427 


2.934 0.662 


0.416 


0.200 


21.4 


4.56 


5752 


4.20 


+0.02 


Mishenina et al ( 2003 ) 














5765 


4.46 


+0.09 


Santos et al J2003H 














5685 


4.26 


+0.07 


Laws & Gonzalez l|2()()l|l 














5760 


4.40 


+0.06 


Deliyannis et al (1200011 














5700 


4.35 


+0.06 


Gonzalez fHMt 














5773 


4.42 


+0.06 


Feltzing & Gustafsson 1)1998^ 














5766 


4.29 


+0.05 


Fuhrmann et al (I1997H 














5752 


- 


- 


di Benedetto, 11998H 














5767 


- 


- 


Alonso et al (1996a) 














5753 


4.33 


+0.06 


Friel et al 1119931 














5753 


4.35 


+0.06 




HD010145 


3.003 0.689 


0.421 


0.212 


36.7 


4.87 


5673 


4.40 


-0.01 


Mishenina et al (2003 ) 














5673 


4.40 


-0.01 





